Shot noise in ferromagnetic single electron tunneling devices 
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Frequency dependent current noise in ferromagnetic double junctions with Coulomb blockade 
is studied theoretically in the limit of sequential tunneling. Two different relaxation processes are 
found in the correlations between spin polarized tunneling currents; low frequency spin fluctuations 
and high frequency charge fluctuations. Spin accumulation in strongly asymmetric junctions is 
£f~) ' shown to lead to a negative differential resistance. We also show that large spin noise activated in 

the range of negative differential resistance gives rise to a significant enhancement of the current 
noise. 
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Recent progress in nanotechnology of magnetic materials renewed interest in spin-polarized electron tunneling in 
magnetic junctions. This interest is additionally stimulated by expected applications of magnetic junctions in field 
sensors, magnetic random access memories (MRAMs), and other microelectronic devices. It is well known since the 
work of JulliereEl, that the junction resistance depends on its magnetic configuration, and this effect, know also as the 
tunnel magnetoresistance (TMR), offers the possibility of magnetic control of the tunneling current. 

A special kind of tunnel junctions are the double barrier junctions with a small central electrode. The increase 
O | in electrostatic energy due to charging this electrode with a single electron can lead to Coulomb blockade of electric 
current below a certain threshold voltage and to Coulomb steps in I-V characteristics at higher voltages. Such 
| junctions, known as single electron transistors (SETs), were extensively studied in the last decade in the limit of 
. nonmagnetic (metallic or semiconducting) electrodes.0 However, it is only very recently when the inteipiay of charge 
t-H ' and spin effects in SETs with ferromagnetic electrodes was studied theoretically and experimentally.lTTj It has been 
shown that the charging effects give rise to periodic modulations of the TMR effect with increasing bias voltage. It 
has been also predicted that cotunneling effects in the Coulomb blockade regime can significantly enhance TMR. 
When intrinsic spin relaxation time on the central electrode (i.e. spin relaxation due to spin-flip scattering within 
the electrode) is much longer than injection time (time between successive tunneling events), then a nonequilibrium 
q\ I spin distribution of electrons is driven by the flowing current. The corresponding spin accumulation can lead to new 
. effects, like negative differential resistance (NDR) or inverse TMR effect. 

A very important parameter, in view of potential applications of the TMR effect, is the noise to signaLratio. 
In the case of nonmagnetic SETs this problem was recently extensively studied in the relevant literature. 0113 It 
has been shown that the noise can provide additional information on the electronic structure, transport properties 
and also on electron-electron interactions. A special issue is the problem of suppression of the zero-frequency shot 
noise by electron correlations. Negative correlations between current pulses, induced by the Pauli principle, can lead 
to complete suppression of the shot noise in quantum point contacts or to reduction of the noise to (l/3)Sp j sson 
(Spoisson — 2el) in metallic diffusive ccmductorsflia Theoretical studies-pf current correlations in the presence of 
Coulomb blockade in SETs by KorotkovEi, Hershfield et al.£H and othersEj showed that the shot noise is reduced to 
{l/2)Sp i SSOn . These predictions were_also confirmed experimentally on quantum point contacts, diffusive conductors, 
' resonant tunneling diodes and SETs.E3~Eil 

However, the situation if ferromagnetic junctions significantly differs from that in nonmagnetic SETs and to our 
best knowledge there is neither experimental nor theoretical work published up to now on shot noise in ferromagnetic 
SETs. As we show in this paper, large spin fluctuations, which occur in ferromagnetic junctions, have a significant 
influence on the shot noise, particularly at low frequencies. We also show that two different relaxation times (spin and 
charge) appear in the frequency dependent current noise. Large spin fluctuations also occur in nonmagnetic junctions, 
but they do not contribute to the current noise. Another important result following from our analysis is the presence 
of an enhanced shot noise (well above 2el) in the NDR range. Similar enhancement was recently observed in resonant 
tunneling diodes. The NDR is due to spin accumulation, which in turn follows from detailed balance of electrons with 
opposite spins tunneling to and off the central electrode. 
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The paper is organized as follows. In Section 2 we describe the formalism which we use for shot noise calculations 
in ferromagnetic SET's with the central electrode small enough so the effects due to charging and discrete energy 
levels are important. The formalism is an extension of the formalism used for nonmagnetic SET's by including spin 
degrees of freedom. In section 3 we analyze the charge-charge and spin-spin correlations as well as the current shot 
noise in ferromagnetic junctions with a nonmagnetic central electrode. The case when one external electrode is made 
of a strong ferromagnet while the other electrodes are nonmagnetic is considered in Section 4. Final conclusions and 
discussion are in section 5. 



II. DESCRIPTION OF MODEL AND METHOD OF CALCULATIONS 



We consider a double junction in which a nonmagnetic metallic grain is separated from two ferromagnetic leads 
(source and sink electrodes) by tunnel barriers with the corresponding spin dependent resistances Ri-\ and on the 
right and i?2| an d R2.1 on the left sides. If the resistances Rj a (j=l,2) are much larger than the quantum resistance 
Rq = h/2e 2 , electronic transpart is dominated by incoherent, sequential tunneling processes and can be described 
within the orthodox approach,E3 in which tunneling is treated in the lowest-order perturbation theory whereas higher- 
order tunneling processes (cotunneling) are neglected. To describe tunneling current and shot noise we extend previous 
descriptions by including the spin degrees of freedom. This generalization is not trivial and leads to effects which are 
absent in the spinless descriptions. 

Electronic transport through the system is govern by the master equation (written in the matrix form) 



dp 



(i) 



for the time evolution of the probability p(Nf, N±;t) to find and A^ excess electrons on the grain. In the limit of 
long intrinsic spin relaxation on the grain (much longer than the time between successive tunneling events) and for 
spin conserving tunneling processes, the matrix M is given by 
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where 



A(N h N l ) = -T+(N h N l )-T- 1 (7V T ,N l )- T+ (iV T , iV ; ) - 1^ (JV t , N{) 
-t+ (tv t , n±) — r- T (7V T ,N l )- T+ (iV T ,N l )- T n (JV T ,N0, 
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(3) 
(4) 
(5) 



The tunneling rates for electrons with spin a, tunneling to (+) and off (— ) the grain through the j-th junction, are 
given byfca 
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where Ef denotes the Fermi energy and the summation runs over discrete energy levels £W of the grain, which are 
assumed to be equally separated with the level spacing AE. Apart from this, 



Vj (N h N l ) = (-iyV^0 + 



eN 
~C 



(7) 



and N = + is the total number of excess electrons on the grain, are the capacitances of the right of and 
left junctions, C is the total capacitance C = C\ + C2, while e and T stand for the electron charge and temperature, 
respectively. The single-electron charging energy is defined as E c — e 2 /2C. When writing Eq.(6) we assumed that 
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electron energy relaxation time is of the order or shorter than the time between successive tunneling events, so 
electrons of a given spin orientation are in thermal equilibrium. 

The electric current in the stationary state can be calculated from the formula 



I = e E [r^ T (iVT,^i)-rr T (iv T ,iv i ) + r+(7v T ,7v i )_r-(iv T ,7v i )]p (7v T ,iv i ) 



= e ]T [r- T (tv t ,N l )- r+ (iv T ,n 1 ) + (jv T , jv ; ) - r+ (tv t , Nl )]/ ( tv t , ^ ) 



(8) 



where the probability p°(N^, N\) is a solution of the equation 6 = Alp . One can also determine the average value of 
any physical quantity X in the stationary state as 



<X>= X^jPiN^Nj 



(9) 



where Xj^. i n 1 is the representation of X in the space of states described by the numbers Nf and N± of excess electrons 
on the grain. _ . 

To analyze fluctuations in the system we extend the generation-recombin&tioti approachE3 for multi-electron channels 
by generalization of the method dejzfiloped for spinlcss electrons in SETlI9~II3. The time correlation function of the 



quantities X and Y is expressed as 



<X(t)Y(0)>= Yl X N , vN[ P(N' v N{;t\N h N l ;0)Y N ^ Nl p°(N h N l ) 

N' v N' l ;N- l ,N l 



(10) 



Here, P(Np iVj; t\N^, N^; 0) is the conditional probability to find the system in the final state with iVj and JV| excess 
electrons at time t, if there was JVj and N± excess electrons in the initial time t=0. This probability is determined 
from the equation 



dP 
~dt 



= MP . 



(11) 



In calculations of the current-current correlation functions one has to include the self-correlation terms as wellEit2l 
According to this procedure the Fourier transform of the charge-charge Snn (upper sign) and spin-spin Smm (lower 
sign) correlation functions are given by 



'jvjv(mm) (w) =4 ]T (N\ ± N{) Re 
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while the current-current correlation function is 

S n (u) = Sff + S c u (") , 
where the Schottky value (the frequency independent part) is 



(12) 



(13) 



qSh _ 



C1C2 



3 - J ' N^N V 

and the second term in Eq.(13) is given by 
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In the present case the Green's function G(u>) = [icol — M] _1 is-dcfini 
in contrast to the previous approaches for the spinless SETsEjEJ, 
the one-dimensional space (i.e. it was tridiagonal) . 



^ed in the two-dimensional space of states (iV-j- , N± ), 
where the corresponding Green's function was in 



III. FLUCTUATIONS AND NOISE IN F-N-F JUNCTIONS 



Let us begin our discussion with a nonmagnetic grain connected to two ferromagnetic leads (F-N-F junction). 
In Fig. la the current-voltage characteristics are presented for the antiparallel and the parallel configurations. We 
assumed there the charging energy E c equal to lO.lmeV and the level spacing AE equal to 3meV. This means that 
for the electronic bandwidth W w lOeV the metallic grain contains N a ~ 10 3 -10 4 atoms (AE ss W/N a ). The 
temperature T = 2.3K is small enough to see the effects due to discrete charging and discrete structure of the energy 
levels (ksT <§; E C ,AE). In Fig. lb the charge and spin accumulated on the grain are shown as a function of the 
bias voltage. The steps in the charge accumulation are related to the steps of the current as they both are caused 
by discrete charging effects. As the charge accumulation is similar in both configurations, the spin accumulation 
only occurs for the antiparallel configuration and is absent in the parallel one. This is because we assumed that the 
ferromagnetic electrodes are made of the same materials, so Ri^/R±i = R21/R21 in the parallel configuration. The 
oscillations of the spin accumulation are due to periodic blockade of the channels for electrons with the spin a =| and 
a =4. With increasing V the mutual balance between the currents through both junctions leads to an increase in the 
spin accumulation and this gives rise to shifts of the chemical potentials for electrons with opposite spins. When V 
increases further, a new channel for spin a ={ opens (in the present case at V « 28mV), which leads to a decrease 
of the spin accumulation as the charge accumulation is almost constant. At V rs 40mV a new channel for electrons 
with a =| opens as well and now both the spin and charge accumulations increase and the next cycle begins. Small 
steps in the curves shown in Fig.l (and also in the following numerical results) is due to discrete structure of the grain 
energy levels. 



In order to study the frequency dependence of the correlation functions Snn(^), Smm{<+>) and Su(uj) we perform 
their spectral decomposition. Accordingly, we express the correlation functions as the sum 



of the components <S*xx(A) m the representation of the eigenvalues A of the matrix M. All eigenvalues A are negative, 
except the eigenvalue A = which corresponds to the stationary solution. The bias dependence of the eigenvalues 
A is presented in Fig. 2 for the spinless and ferromagnetic SETs (the upper and lower parts, respectively). In the 
latter case the structure is much more complex than in the former one because the spin degeneracy is removed in 
ferromagnetic SETs. There is also hybridization of the eigenvalues, especially strong in the low frequency region and 
for voltages close to the steps in the I-V curves. 

Calculations of the noise were performed numerically for the space of states (JVf, iVjJ dynamically increasing with 
V, and the convergence of the results was checked for the large space of states, which in an extreme situation included 
247 states. The results of the spectral decomposition were compared with those obtained_hy direct calculation of the 
Green's function G(u>) for some small values of u> (by inversion of the matrix [iuil — M])Ji3 

An example of the spectral decomposition of the charge-charge, Snn(^)i and spin-spin, Smm(w), correlation 
functions is shown in Fig. 3. The amplitude of the spin noise <S*mm(A) is two orders of magnitude larger than the 
amplitude of the charge noise S , jvjv(A). From Fig. 3 follows that only a few components are significant; the others 
are exponentially small and can be omitted. It would be rather difficult to extract experimentally the contributions 



A. Charge and spin fluctuations 




A<0 



(16) 
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corresponding to A's, which are close to each other. In such a case one can approximate Snn(u) with a single 
Lorentzian curve by introducing an effective relaxation time, 

S NN (U) = 4 7~charge 

1 + ^charge 

as it would be for poissonian processes of independent tunneling events. Here, var(iV) =< iV 2 > — < N > 2 is the 
variance of N. In Fig. 4 we present the characteristics for the charge noise: the variance var(TV) and the inverse of the 
effective relaxation time l/r c ^ Qrge . The voltage dependence of 1/T c h arge is a periodic function, which has minima at 
the steps of the charge accumulation and maxima at the points corresponding to the minima of var(iV). The value 
of 1 / T c harge is an order of magnitude larger than the average frequency of electron transfer from the source to sink 
electrodes, l/r tr = I / e w 1GHz. This is because the current / is determined by the more resistive junction, for which 
the electron tunneling frequency is close to 1/T tr , while 1/T c ^ arge (~ E c /mxa[Rj a )) is determined by the frequency of 
electron tunneling through the less resistive barrier. 

Since we assumed the level spacing AE < E c , we can expect large spin fluctuations (var(M) larger than var(TV)) 
and a long spin relaxation time (r sp i n > T cna rge)- The analysis of the spin noise has been performed in a similar way as 
of the charge noise. The corresponding results presented in Fig. 5 show that in the case considered l/r sp i„ is 20 times 
smaller than 1/ Tcharge- Since l/r sp i„ ~ Ai?/max(i? :((T ), its value can be even smaller for a larger number of atoms 
in the grain (N a ss W/AE). There are also large differences in var(M) and l/r sp i n in the antiparallel and parallel 
configurations. The spin fluctuations are significantly smaller in the antiparallel configuration than in the parallel 
one. This difference is due to large spin accumulation in the antiparallel configuration, which shrinks the space of 
available states for the spin fluctuations. The small steps in var(M) and 1/T sp i„, clearly seen in Fig. 5, particularly 
in the antiparallel configuration, result from the step-like changes of the spin accumulation on the grain. Using the 
spectral decomposition procedure we found that the dominant contribution to the spin noise comes from the largest 
eigenvalue Ai. This can be concluded from Fig. 5b, where the exact value of l/r sp i„ (thick solid and dashed curves) is 
compared with the contribution corresponding to Ai (thin solid and dashed curves). 

B. Current noise 

The current shot noise has been calculated from the formulas (|ll|)-(|l5|). Figure 6a shows the bias dependence of the 
zero- frequency current noise Su(lu = 0). The current noise is always smaller in the antiparallel configuration than in 
the parallel one. This is because in the presence of spin accumulation the amplitude of fluctuations is smaller for the 
antiparallel configuration. In Fig. 6b Su(oj = 0) is split into two parts; a frequency independent component Sfj 1 and a 
contribution Sjj(uj = 0) arising from the frequency dependent part of the current noise (see Eq.(13)). The component 
Sfi is almost constant sa 2eI(C\ + Cf )/C 2 at the plateaux of the I — V curve and increases with opening of new 
channels. Dynamical correlations between the currents are described by Sjj(lu). Its value in the limit uj — > can be 
positive between the I — V steps and negative when new channels become open. This is evident for the antiparallel 
configuration at V ~ 26mV when opening a tunneling channel for electrons with u =[ leads to negative dynamical 
correlations. This effect is almost compensated by an increase in Sfj and therefore, one gets only a small reduction 
of the current noise Su{uj = 0). 

A large peak of Sfp seen in Fig. 6b at the threshold of the Coulomb blockade (V ss 14mV) occurs in any asymmetrical 
SET device. It is caused by a rapid increase of a stream of tunneling electrons onto the grain, which appears when 
the conducting channels open. Since the outgoing channel has larger resistance the charge accumulation increases and 
part of electrons is pushed back to the source electrode. Frequent tunneling of electrons to and from the grain leads 
to an enhancement of the charge noise as well as of the current noise (see also Fig. 4). The reduction of Su(uj = 0) at 
V ~ 14mV and V ~ 40mV is the effect of the opening of new current channels at the steps of the I — V curves. 

Frequency dependence of the Sjj(uj) component of the current noise is presented in Fig. 7a for different voltages. 
One can clearly distinguish between two different relaxation processes in the current noise; respectively in the low 
and high frequency ranges. An interesting case is that for V = 36 mV (and also V = 42 mV), where anticorrelations 
between currents occur in a low frequency limit (with Sfj < 0) and correlations for higher frequencies (Sjj > 0). 

In order to determine the effective relaxation time tjj we have performed a detail analysis of Sjj(uj) as a function 
of uj. The shape of Sjj{u)) suggests that there are only a few relaxation processes with well separated relaxation times 
tjj. Therefore, one can assume that Sjj(uj) is a superposition of a few Lorentzian curves. The relaxation time for 
each process is given by 1/ Tjj = 3ojq , where ujq corresponds to the maximal slope of Sjj and can be determined from 
cPSfj/duj 2 — 0. The results are presented in Fig. 7b by open circles. One can distinguish two distinct relaxation times, 
one in the high and another one in the low frequency ranges. In a wide voltage range the corresponding relaxation 
times are very close to the effective relaxation times for the charge and spin noise (compare with the solid curves in 
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Fig. 7b). Some deviations occur for voltages where new tunneling channels become open. In a certain voltage range 
there are three relaxation times (e.g. at V ~ 36mV - see also the curve with squares in Fig. 7a). This analysis shows 
that both charge and spin fluctuations are relevant for the shot noise of the current in ferromagnetic SETs. 

C. Paramagnetic limit 

An interesting limit of the F-N-F junction discussed above is the case where all electrodes are nonmagnetic. It is 
important to note that large spin fluctuations are present in such the case as well. Fig. 8 shows characteristic frequency 
dependence of Sjj as well as of Snn and Smm- Although the amplitude of the spin noise is very large (in this case it 
is four orders of magnitude larger than the charge noise), it is not seen in the current noise. In a nonmagnetic system 
the tunnel resistances for both spin directions are the same. This implies that the components of the current-current 
correlation functions obey the relation Sf j. (A) = S/ j. (A). For those eigenvalues A which correspond to the spin 
noise we have also Sf j (X) = —Sfj(X) < 0. This means that the tunneling events with the frequency 1/T sp i n 
are anticorrelated for electrons with the same spin and correlated for electrons with the opposite spin. Therefore, the 
spin component of the total current noise is completely compensated for each junction in a nonmagnetic SET. 

The asymmetry between the tunneling channels for electrons with the opposite spins leads to activation of the spin 
component in the current noise. We extracted from Sjj the components Sjj charge and Sjj spin corresponding to the 
charge and the spin noise, respectively. Their values in the zero-frequency limit are determined from Sjj charge (u> — 
0) = -4^ A<Ao Su(X)/X and Sjj spin (u = 0) = -iJ2x <X<o S nW/\ where A o = ^r charge r spm . The results 
are presented in Fig. 9 as a function of the parameter p, where p\ = pi with pj — {Rj-\ — Rji) I {Rtf + Rji) f° r J = 1, 
2. It can be seen that the charge component is almost constant whereas the spin component increases with p and for 
p — > 1 can be much larger than the charge component. 

IV. ENHANCEMENT OF THE NOISE IN F-N-N JUNCTIONS 

Let us consider now the electron transport through a nonmagnetic metallic grain connected to a nonmagnetic lead 
on one side and to a ferromagnetic lead on the other side (F-N-N junction). When the ferromagnetic electrode is 
a strong ferromagnet with only one spin subband (say for spin a =1) partially filled, then the electrons with a =f 
can tunnel onto the grain from the nonmagnetic electrode, but the outgoing a —\ channel from the grain to the 
ferromagnetic electrode is blocked (Ri\ — oo). Such a double junction can be considered as a system of coupled 
SET (for a ={) and a single electron box (SEB) (for a —1). Nonmagnetic coupled SET-SEB devices were recently 
analyzed-,ti«3oretically aaad it was shown that if the coupling is strong, one can have a negative differential resistance 
(NDR).E!lc2l Heij et alxl constructed recently the SET-SEB device based on two electrostatically coupled metallic 
grains, and found NDR indeed. In our case we have only a single grain, but two coupled electron channels. In the I-V 
curve shown in Fig. 10 we see NDR with a sequential drop of the current. With increasing V the number of electrons 
with a =| on the grain increases (see the spin accumulation in Fig. 10b), which results in an increase of the local 
chemical potential for a =\ electrons and a drop of the current. The maximum drop of I occurs at the maximum of 
the spin accumulation (at V ~ 26mV). Opening of the next charge channel (N — > N + 1) for conducting electrons 
(a =1) reduces the spin accumulation and increases the current. Small steps in the charge and spin accumulations 
are the effect of opening of spin channels, i.e. for (iV-j-, JVj.) — > (iVj ± 1, iVj_ =F 1). 

The noise in magnetic SET-SEB devices has been analyzed with the same method as before and the results for the 
charge and spin noise are presented in Fig. 11. The magnetic fluctuations increase when a new channel opens and then 
they drop again with increasing voltage. The peaks in var(Af) correspond to the maxima of the effective relaxation 
time T S pi n (the minima of 1/T sp i n shown in Fig. lib). The spin noise has a significant influence on the current noise 
shown in Fig. 12. The large increase of the current noise (much over the value Sp i SSOn ) takes place for the range of 
NDR, which occurs for voltages corresponding to large spin fluctuations. We performed the spectral decomposition 
of the current noise Sjj(tu). The voltage dependence of the eigenvalues A did not show any peculiarities in the range 
of NDR. The spectral weight is, however, shifted to lower frequencies. The frequency dependence of Sjj(u>), shown 
in Fig. 13, exhibits strong increase of the low frequency contribution for the voltages where NDR occurs. Outside the 
NDR range the low frequency-component of the current noise is minor. 

Recently, Iannaccone et al.E3 studied resonant tunneling through a double barrier GaAs structure and showed an 
enhancement of the shot noise in the region o£-NDR. The NDR effect was due to charge accumulation in the central 
well (the effect predicted by Ricco and AzbeEj) . The, enhancement of the shot noise was also observed in resonant 
tunneling through a double barrier in magnetic field^J although NDR_was a result of changes in electronic structure 
in high magnetic field (formation of Landau levels). Iannaccone et al.E3 analyzed the enhancement of the shot noise 
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in the framework of the generation-recombination approacht§c3 c3 in terms of the effective generation (g) and the 
recombination (r) rates-aiidr-their relaxation times r g and r r , respectively. Such the approach is equivalent to two 
channel approximation.llH~t3E3 In the present studies we used the method which takes into account multichannel (and 
multilevel) tunneling processes and therefore we have insight into the microscopic processes related to the tunneling 
current and its noise. 



V. FINAL REMARKS 



In summary, we have performed theoretical analysis of the frequency dependent current shot noise in ferromagnetic 
tunnel junctions in the limit of sequential tunneling. We showed that apart from charge fluctuations there are also 
strong spin fluctuations, which influence the current noise. The spin noise is significant at low frequencies, while the 
charge noise at higher frequencies. We predict that two distinct relaxation times r sp i n and T c / large should be seen in 
frequency measurements of the current noise in magnetic tunnel junctions. 

We have also studied the magnetic SET-SEB device. The spin accumulation effect can lead in this case to a negative 
differential resistance, which occurs only when the spin accumulation increases, i.e., when the chemical potential of 
the conducting electrons increases due to accumulation. Spin fluctuations are then activated, which leads to a strong 
enhancement of the low frequency component of the current noise. 

The numerical results have been calculated at a low temperature (kpT -C E C ,AE), where all features due to 
discreteness of the system are clearly seen. We have also performed calculations for higher temperatures. The small 
steps, corresponding to the opening of the magnetic channels, are washed out when T — * AE/kp. However, the charge 
and spin fluctuations, var(iV) and var(M), and the corresponding relaxation times, r c ^ arge and T sp i n , remain almost 
unchanged. In F-N-F junctions (as in Fig. I) the tunnel magnetoresistance, defined as TMR = Ip/Iap ~ I (where Ip 
and Iap are the currents in the parallel and in the antiparallel configurations), decreases with increasing temperature 
(see also Ref.[ ^|). For T > AE/kp, the physical quantities become less dependent on the voltage. The values of 
var(TV) and var(M) increase with T, as expected. The charge and spin fluctuations are still well separated in the 
different frequency ranges. At T = 34. 8K the zero-frequency current noise Sjj(lu = 0) ~ 0-8Sp O i SSOn . The reduction 
is due to Coulomb correlations, which are still large as T < E c /kp — II7K. The component Sfj is, however, much 
larger than Sp i SSOn , whereas the frequency dependent part Sjj(u>) is negative with a large amplitude. Its low and 
high frequency components, corresponding to the spin and the charge fluctuations, are also well defined. 
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Voltage [mV] 



FIG. 1. Voltage dependences of the current (a) and charge (solid and dashed curves) and spin (curves with circles and 
crosses) accumulation (b) in the antiparallel and the parallel configurations. The parameters assumed here are: i?2T — 240MQ, 
R 2i = 60Mfi, Rn = 2MQ. and R n = 8MQ for the antiparallel configuration (fl 2T = 240Mft, R 2i = 60Mf7, R^ = 8MD. and 
R n = 2MQ for the parallel configuration), C 2 = 6.6aF, Ci = 1.32aF, AE = 3meV and T = 2.3K. The inset in part (a) shows 
the scheme of the system (electrons flow from right to left). 
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1 20 30 40 50 60 
Voltage [mV] 

FIG. 2. Voltage dependence of the eigenvalues A of the matrix M for the spinless (the upper part) and ferromagnetic (bottom 
part) SETs. The parameters assumed in (a) are: R 2 = 240Mf2, Ri = 2MQ, C 2 = 6.6aF, Ci = 1.32aF and AE = 3meV, while 
in (b) are the same as in Fig.l for the antiparallel configuration. The space of states is confined to — 5 < N < 5 for the spinless 
SET and -1 < N < 4, -8 < M < 8 for the ferromagnetic SET. 
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FIG. 3. The components of the spectral decomposition of the charge-charge and the spin-spin correlation functions: S'jvjv(A) 
(dots) and Smm(A) (triangles), determined for V = 25mV. The other parameters are as in Fig.l for the antiparallel configura- 
tion. 
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FIG. 4. Voltage dependence of the variance of charge fluctuations, var(iV), (a) and of the inverse effective relaxation time, 
1/T c harge, (b) for the parallel and antiparallel configurations. The parameters are the same as in Fig.l. 
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FIG. 5. Voltage dependence of the variance var(M) of the spin fluctuations (a) and of the inverse effective relaxation time 
1/Tspin (b) in the antiparallel and the parallel configurations. The parameters are the same as in Fig.f. The thin (solid and 
dashed) curves in part (b) present contribution corresponding to Ai , where Ai is the largest eigenvalue of the matrix M. 
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FIG. 6. Voltage dependence of the current shot noise at u> = (a) for the system defined in Fig.l. In part (b) Sn(ui = 0) is 
split into two components: Sff (upper curves) and S c u {uj = 0) (lower curves). 
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FIG. 7. (a) Frequency dependence of the current shot noise £//(<*;) in the system defined in Fig.l for the antiparallel 
configuration, (b) Inverse of the relaxation time for the current noise (open circles) as a function of V. For comparison the 
inverse of the effective relaxation times l/r c h ar9<3 and 1/T sp i„ for the charge and the spin noise are shown by the solid upper 
and the lower curve, respectively (the same curves as the solid ones in Fig. 4b and Fig. 5b). 
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FIG. 8. Frequency dependence of the current noise Sij(u>) (a), and of the spin and charge noise (b) at V = 26mV in a 
nonmagnetic junction. The parameters assumed are: = Rii — 2MQ, R21 = R21 = 60Mf2, C2 = 6.6aF, C\ = 1.32aF, 
AE = 3meV and T = 2.3K. 
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FIG. 9. Two components of «S/j(w = 0) corresponding to the charge and spin noise plotted as a function of 
p = (Rii — Rii)/(Rii + Ri.i) = (i?2T — #2|)/(-R2T + R2.i)- The other parameters assumed here are: R21 = 60Mf2, 
i? 2T = (1 -P)R*1, Ril = 2Mfi, fli T = (l+p)/{l-p)R n , Gi = 6.6aF, Ci = 1.32aF, AE = 3meV and T = 2.3K. 
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FIG. 10. The I — V curve (a) and the charge and spin accumulation (b) for a F-N-N junction. The parameters assumed are: 
i?2t = oo, R21 = 9Mfi, Rn = R n = 2MQ, C 2 = 6.6aF, Ci = 1.32aF, AE = 3meV and T = 2.3K. The inset shows the scheme 
of the system (electrons flow from the nonmagnetic to the ferromagnetic electrode). 
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FIG. 11. The variance of the charge and spin fluctuations (a) and the inverse effective relaxation times l/r c harge (solid curve) 
and 1 / Tspin (dashed curve) (b) as a function of the bias voltage for the system as in Fig. 10. 
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FIG. 12. Voltage dependence of the zero- frequency current shot noise (a) in the system defined in Fig. 10. In part (b) two 
components of the shot noise Sfj 1 (solid curve) and Sjj(oj = 0) (dashed curve) are presented. 
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FIG. 13. Frequency dependence of the current noise Sfj(u)) in the system as in Fig.10 for indicated voltages. The curves for 
V — 21.4mV and V = 25.1mV correspond to the maxima of the current shot noise seen in Fig. 12 in the range of NDR. 
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